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We show through simulation that quantum interference in non-sequential double ionization can 
be used to control the recollision with subsequent ionization (RESI) mechanism. This includes the 
shape, localization and symmetry of RESI electron-momentum distributions, which may be shifted 
from a correlated to an anti-correlated distribution or vice versa, far below the direct ionization 
threshold intensity. As a testing ground, we reproduce recent experimental results by employing 
specific coherent superpositions of excitation channels. We examine two types of interference, from 
electron indistinguishability and intra-cycle events, and from different excitation channels. 


Correlation and anti-correlation have been extensively 
studied in strong-field, laser-induced nonsequential dou¬ 
ble ionization (NSDI). In particular with the Cold Tar¬ 
get Recoil Ion Momentum Spectroscopy (COLTRIMS) 
technique, information about the electron momenta has 
become experimentally accessible since the early 2000s. 
Several features in these distributions provide informa¬ 
tion about the type of interaction by which NSDI occur, 
and the physical mechanisms behind it. It is commonly 
accepted that NSDI results from the laser-induced in¬ 
elastic recollision of an electron with its parent ion [1]. 
If the driving-field intensity is high enough, upon recolli¬ 
sion the first active electron releases a second electron by 
electron-impact (El) ionization. In contrast, for the so- 
called “below-threshold” intensities, the kinetic energy 
transferred from the first electron to the core is not suffi¬ 
cient to free the second electron. Instead, the recolliding 
electron imparts only enough energy so that the second 
electron is excited, and then ionized with a time delay. 
This mechanism is known as recollision-exctiation with 
subsequent ionization (RESI). 

Throughout the years, the prevalent view is that, in 
El, the two electrons will exhibit correlated momenta, as 
a consequence of their being released simultaneously. In 
contrast, for RESI, back-to-back emission will occur due 
to the time delay between recollision of the first electron 
and ionization of the second electron. This view has been 
backed using classical models, which have reproduced 
many of the key features encountered in experiments (for 
reviews see [2]). Such models also exhibit excellent agree¬ 
ment with the outcome of other methods, such as the 
full solution of the time-dependent Schrodinger equation 
[3], or the strong-field approximation (SFA) [4-6]. The 
above-mentioned studies, performed for El, suggested 
that quantum mechanical features such as interference 
will not survive integration over momentum components 
perpendicular to the laser-field polarization, which is the 
typical scenario in NSDI experiments. These conclusions 


were then extrapolated to RESI without much evidence. 

Classical models have also successfully reproduced the 
anti-correlated behavior observed in early RESI experi¬ 
ments [7-10]. However, recent RESI experiments have 
found that electron-momentum distributions may occur 
in a variety of shapes, from cross-shaped distributions lo¬ 
calized along the axis [11] to probability densities spread 
between all four quadrants of the plane spanned by the 
parallel electron momentum components [12, 13]. These 
findings contradict the simple explanation that a time 
delay will lead to anti-correlated electron momenta. 

A wealth of shapes has also been observed in many the¬ 
oretical studies as diverse as the SFA and similar meth¬ 
ods [14-16] to classical trajectory computations [17-19]. 
If quantum interference is absent, the SEA retains four¬ 
fold symmetry for RESI distributions, in agreement with 
rigorous momentum constrains [15]. Recently, however, 
SEA computations have shown that this symmetry can 
be broken, if the interference stemming from different 
excitation channels is incorporated [20]. Therein, it has 
been argued that inter-channel interference is paramount 
for obtaining anti-correlated RESI distributions. 

In this Letter, we show that one may obtain corre¬ 
lated, anti-correlated, cross- or ring-shaped RESI distri¬ 
butions in an SEA computation, by choosing appropriate 
coherent superpositions of excitation channels. We re¬ 
produce recent experimental results in which dramatic 
variations in the shapes of RESI distributions have been 
observed [12]. The features encountered are related to 
the interplay between two types of interference, involving 
(a) events which are displaced by half a cycle and those 
present due to the symmetry of indistinguishable elec¬ 
trons; and (b) different channels of excitation for the sec¬ 
ond electron. These effects individually break the four¬ 
fold symmetry of the momentum distributions, and may 
be used to manipulate the electron-electron correlation. 

In the SFA and in atomic units, the RESI transition 
amplitude related to the c-th excitation channel reads 
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describes the process in which an electron, initially at 
a bound state of energy —E\ g , leaves at t", returns to 
the core at t' with intermediate momentum k and ex¬ 
cites a second electron from a state with energy —E 2g to 
a state with energy —E^*. The first electron is released 
at t' with momentum pi, while the second electron is 
freed at a subsequent time t with final momentum p 2 . 
The prefactors , V^], kf; and Vpf e are related to the 
ionization of the first electron, the recollision-excitation 
process and the tunnel ionization of the second electron, 
respectively. They contain all information about the in¬ 
teractions, which are chosen as long range, and electronic 
bound states [14, 15]. We compute the transition ampli¬ 
tude (1) using the steepest descent method, as described 
in [21]. We use the notation p n || and p„j_, n = 1,2, for 
the momentum components parallel and perpendicular 
to the laser-field polarization, respectively. We employ 
a monochromatic field, which is a good approximation 
for longer pulses, and reasonable for short pulses if the 
carrier-envelope phase is averaged out. 

Due to momentum constraints attributed to the rescat¬ 
tering of the first electron and ionization of the sec¬ 
ond electron occurring most probably near field cross¬ 
ings and crests, respectively, the distributions will be 
located around (pi||,_P 2 ||) = (±2\/0p,0). Considering 
the symmetry of displacement by half a cycle and par¬ 
ticle exchange of the system leads to four “events” lo¬ 
cated around (— 2y/Up, 0), (2y/Up, 0), (0,2 yJTJp) and 

(0,-2 y/Up). We label these event amplitudes as M^ c \ 

Mr°\ Mu ^ and M^ c \ given they occupy the left, up¬ 
per, right and lower regions of the parallel momentum 
plane respectively. Quantum interference of events oc¬ 
curs mainly in the overlap of the above regions, which is 
located around the diagonals pi\\ = ±p 2 ||- 

For a single channel, the coherent sum of such ampli¬ 
tudes, integrated over the transverse momentum compo¬ 
nents reads 


tt {c) (pi\\,p 2 \\) 


d 2 pi±d 2 p 2 ± 


M (c 


( 3 ) 


with M = M^ + Mr^ + M^ If the events 

are summed incoherently, the amplitudes are replaced by 
probabilities, i.e., \M^ C ' I \ 2 + |Mr C) | 2 + |M„ C ^| 2 + \M^\ 2 is 
employed instead. If more than one channel is involved, 
Eq. (3) is generalized to 


tt(pi\\,P2\\) = / d pi±d p 2 ± 
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Channel Interference 



FIG. 1: (Color online) Schematic representation of the event 
and inter-channel interference for RESI (panels (a) and (b), 
respectively). The shaded regions in panel (a) represent re¬ 
gions where M[ c \ Mr C \ Mu'* and M^ c * are substantial, and 
the arrows in panel (b) indicate the different excitation path¬ 
ways for the second electron, according to Table I. 


Channel 

Transition 

E 2e (a.u.) 

1 

3s3p 6 (3s ->• 3p ) 

0.52 

2 

3p 5 3d (3 p ->• 3d) 

0.41 

3 

3p 5 4s (3p —> 4s) 

0.4 

4 

3p 5 4p (3p —> 4 p) 

0.31 

5 

3p 5 4d (3p —> 4d) 

0.18 

6 

3p s 5s (3p —> 5s) 

0.19 


TABLE I: Excitation channels for Ar + employed in this 
work.From left to right, the columns give the number asso¬ 
ciated with the channel, the electronic configurations for the 
sub-levels involved in the excitation and the absolute value 
I? 2 e of the excited-state energy, respectively. The excitation 
pathway for the second electron is given in brackets. 

Eq. (4) assumes that each excitation represents a path, 
which the second electron can take. Hence, the am¬ 
plitudes corresponding to each path must be summed. 
One may include different amplitudes or phases for each 
channel, which would model channel selection or ac¬ 
count for phase effects. This yields a more general sum 
/ d 2 pi±d 2 p 2 ±\J2 c N c e’ lcl,c M^ c *\ 2 . In Fig. 1, we provide a 
schematic representation of the event and inter-channel 
interference, for the excitation channels in Table I. 

Experimental results from [12] of RESI in argon show 
that increasing pulse lengths from few to many cycle 
pulses causes a transition from cross-shaped to slightly 
anti-correlated, correlated or ring-shaped distributions. 
Our results, displayed in Fig. 2, exhibit many of these 
features. Specifically, we can model the transition from a 
cross to an anti-correlated distribution, as shown in the 
top row of Fig. 2. Moreover, we can see in panels e) 
and f) that the distributions go from anti-correlated to 
correlated for increasing laser-field intensity. 

These effects are achieved by using coherent superpo¬ 
sitions of the excitation channels given in Table I, using 
the relative weights and phases in Table. II. The phases 
have been chosen to achieve an anti-correlated pattern in 
panels c). If a pattern is optimized to be anti-correlated 
at one driving field strength, changing this will cause it 
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FIG. 2: (Color online) RESI distributions for argon (E\ g = 
0.58, E 2 g = 1.02 a.u.) computed using different coherent su¬ 
perpositions of the excitation channels in Table I. The phases 
and weights employed in these superpositions are provided in 
Table II, and have been chosen in such a way as to reproduce 
the experimental data from [12]. In the first row, panels a)- 
d), the intensity is 1 x 10 14 W/cm 2 with to = 0.057 a.u. and 
a ponderomotive energy of U p = 0.22 a.u. The intensities in 
the second row are as follows: e) I = 0.8 x 10 14 W/cm 2 with 
Up = 0.18 a.u., f) and g) / = 0.8xl0 14 W/cm 2 with U p = 0.26 
a.u. and h) I = 1.4 x 10 14 W/cm 2 with U p = 0.27 a.u. The 
dashed lines in the figure indicate the diagonals = ±_p 2 || - 
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TABLE II: Coherent superpositions employed in Fig. 2. The 
letters in the first column indicate the panels in Fig. 2 for 
which a specific superposition have been employed. From 
the second to seventh column, the numbers in the first row 
indicate the excitation channel in Table I, and the excited 
state of the second electron is given in brackets. The numbers 
N c e l<t>c give the weight and the relative phase for each channel. 

to flip to be correlated. Hence, in panels f)-h) we obtain 
correlated patterns, in agreement with [12]. One should 
note that, despite this myriad of shapes, all intensities 
employed in this experiment are well within the below- 
threshold regime, as the maximal kinetic energy of the 
returning electron, 3.17 U p , is significantly lower than the 
second ionization potential E 2g [25]. Hence, RESI is the 
prevalent NSDI mechanism. Similar superpositions may 
be used to reproduce the results in [7, 9, 10]. 

The features observed in Figs. 2(a) to (d) mark a 
change from a regime in which excitation to s states is 
prevalent, to a scenario in which a coherent superposi¬ 
tion of p and d states dominates. This statement can be 
inferred from Fig. 3, which shows very different shapes 
for different channels of excitation. These differences 
stem from the prefactor V P2e related to the ionization 
of the second electron. For p and d states, V P2e has an¬ 
gular nodes, which prevent the distributions from being 
located along the axes. In contrast, for s states only ra¬ 
dial nodes are present, so that the electron-momentum 


distributions will be cross-shaped [14, 22]. Furthermore, 
none of the coherent distributions in Figs. 2 and 3 exhibit 
the fourfold symmetry obtained if incoherent sums are 
used. This symmetry breaking occurs already for a single 
channel, with fringes along and parallel to the two diag¬ 
onals. These fringes may be understood if one consid¬ 
ers a single-channel coherent sum (3) and its incoherent 
counterpart, neglecting the prefactors (see Fig. 4). They 
come from the interference of the amplitudes associated 
with electron indistinguishability. The diagonal maxima 
and minima satisfy the condition p 1 y = p 2 \\ ± |<5|, with 
|(5| ~ ojnir/(2y/Up) : where n is even or odd, respectively. 
The anti-diagonal fringes are four times narrower. These 
fringes are obtained analytically by maximizing the inte¬ 
grand of Eq. (3), and have been derived elsewhere [22]. 

In summary, we have shown that two types of quan¬ 
tum interference play an important role in reproducing 
results similar to those in experiments. Event interfer¬ 
ence will break the four-fold symmetry previously seen 
in the SFA, even when only considering a single chan¬ 
nel of excitation. Inter-channel interference can be used 
to maximize anti-correlation, correlation, and alter the 
shape of the electron-momentum distributions, and thus 
create all distributions found experimentally [7, 9, 10, 12]. 

This brings additional insight into many experimental 
studies [7, 9, 12], in which correlated distributions found 
for increasing driving-field intensity have been attributed 
to direct ionization, despite still being below the thresh¬ 
old intensity [10]. Our results suggest that this could in 
fact be RESI, provided an appropriate coherent superpo¬ 
sition of channels is chosen. Furthermore, recent experi¬ 
ments have shown by using a two-color field and changing 
the relative phase of the colors the momentum distribu¬ 
tions for NSDI below the electron impact threshold can 
be manipulated from being anti-correlated to correlated 
[23]. The phases of each channel and event have a strong 
dependence on the field as they derive from the action. 
Hence, it would not be unreasonable to assume this effect 
came from the interference we have described. 

However, one may ask why very short pulses favor 
excitation to s states, while longer pulses favor p or d 
channels. We have found that the relative s-state con¬ 
tributions decrease with the driving-field intensity, while 
those associated to the remaining channels remain rela¬ 
tively stable (not shown). For a few-cycle pulse, there 
are few dominant events near the peak of the pulse and 
other events near smaller crests. In contrast, the inten¬ 
sity across a long pulse will be more uniform. This may 
reduce the s contributions for the same peak intensity. 

Finally, given the important role quantum interference 
has in RESI, one should consider the implications this 
has for classical models. We have shown that interfer¬ 
ence hugely increases the SFA’s predictive power and 
can reproduce many of the features seen in experiment. 
Thus, classical-trajectories should be used with caution. 
Nonetheless, a large density of states may lead to a quasi¬ 
continuum, which could give rise to quasi-classical wave 
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FIG. 3: (Color online) One-channel RESI electron momentum distributions computed for argon with all prefactors included, 
for the parameters in Table I. The laser intensity and frequency are I = 4.56 x 10 13 W/cm 2 and u> = 0.057 a.u., respectively. 
Each panel has been labeled with the corresponding channel number (top left), as stated in Table I, and the excited state of the 
second electron (top right). The top and bottom rows display a coherent and incoherent sum of events, respectively. Angular 
and radial nodes from V P2e have been marked on the bottom row with red lines and green circles, respectively. 



FIG. 4: (Color online) One-channel coherent and incoherent 
sums of all events, computed neglecting the prefactors 
for the same field parameters as in the previous figure [panel 
a) and b) respectively], and the RESI channel 1 in Table I. 
The diagonals are indicated with the red lines in the figure. 

packets where interference would play less of a role. Our 
implementation of the SFA currently neglects broaden¬ 


ing of states caused by the field. However, recent stud¬ 
ies of the RESI phase-space dynamics have revealed a 
highly confined region that can be associated with trap¬ 
ping in an excited state [24]. This would justify using 
discrete bound states, and would render interference im¬ 
portant. This provides a large scope for RESI to be used 
as an attosecond-imaging tool in order to probe and re¬ 
construct the intermediate state of the second electron. 
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